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SUMMARY 


A  small  deflection  plate  theory  is  developed  which  includes  higher  order 
approximations  to  the  effects  of  transverse  shear.  The  governing  equa¬ 
tions  and  associated  boundary  conditions  are  developed  by  a  variational 
method.  The  theory  is  sufficiently  general  to  encompass  anisotropic, 
laminated,  composite,  circularly  curved,  cylindrical  plate  elements  of 
nonsymmetric  cross  section.  Specific  applications  are  made  to  the 
bending  and  the  axial  compression  buckling  of  simply  supported  flat 
plates  and  to  the  axial  compression  buckling  of  curved  plates  and 
shells  having  classical,  simple  supports.  It  is  shown  that  the  effects 
of  interlaminar  shear  can  be  significant  in  composite  constructions 
which  use  very  high  strength  fiber  reinforcement. 
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FOREWORD 


The  work  reported  herein  constitutes  a  portion  of  a  continuing  effort 
being  undertaken  at  Stanford  University  for  the  Eustis  Directorate, 
U.3.  Army  Air  Mobility  Research  and  Development  Laboratory  under 
Contract  DAAJ02-70-C-0075  (Project  lF06ll02A33F)  to  establish  accurate 
theoretical  prediction  capability  for  the  static  and  dynamic  behavior 
of  aircraft  structural  components  using  both  conventional  and  uncon¬ 
ventional  materials.  Predecessor  contracts  supported  investigations 
which  3ed,  in  part,  to  the  results  presented  in  References  10,  11,  lk, 
and  15. 
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INTRODUCTION 


The  desire  for  significantly  increased  efficiency  in  aerospace  structures 

has  encouraged  the  development  of  high-strength  fiber  composites  for  use 

in  laminated  plate  and  shell  elements.  However,  to  realize  the  potential 

inherent  in  composite  construction,  the  theoretical  prediction  capability 

must  be  improved  in  order  to  match  laboratory  experience.  Plate  and  shell 

theories  which  deal  selectively  with  only  some  of  the  characteristics  of 

laminated,  composite  constructions  are  becoming  inadequate  for  the  very 

high  strength  fiber  composites  being  developed.  In  the  case  of  laminated 

1-14- 

plates  and  shells,  researchers  have  been  concerned  primarily  with  the 
anisotropic  nature  of  layered  constructions  under  the  assumption  that  the 
interlaminar  shear  effects  are  negligible.  Whereas  this  assumption  is 
reasonably  valid  for  composites  with  relatively  high  ratios  of  interlaminar 
shear  to  laminae  extensional  moduli  (for  example,  glass-epoxy),  there  is 
reason  to  question  this  assumption  for  seme  of  the  newer  composites  (for 
example,  boron-epoxy)  wherein  relatively  low  ratios  of  these  moduli  exist. 

The  literature  on  the  bending  and  buckling  of  laminated,  composite  plates 
and  shells  may  be  considered  to  begin  with  the  work  of  Smith  (Reference  l). 
Therein,  Smith  considered  the  bending  of  a  two-layered  rectangular  plate 
with  orthotropic  layers  oriented  so  that  the  natural  axes  of  each  layer 
made  equal  but  opposite  angles  with  the  coordinate  axes.  Plates  and  shells 
of  this  construction  are  referred  to  as  being  of  the  Smith  type.  By  over¬ 
looking  coupling  between  bending  and  membrane  stresses,  Smith  concluded 

erroneously  that  a  plate  of  this  construction  behaves  in  the  same  manner 

2 

as  the  conventional  orthotropic  plate.  Reissner  and  Stavsky  formulated 
the  composite  plate  problem  in  terms  of  an  Airy  stress  function  and  the 
lateral  deflection.  For  the  Smith-type  plate,  they  established  the  ex¬ 
istence  of  interaction  between  bending  and  extensional  stresses  in  plates 
under  direct  and  lateral  loads,  respectively.  In  laminated  plates  sym¬ 
metrical  about  a  median  surface,  this  coupling  vanishes.  In  addition, 
Reissner  and  Stavsky  noted  from  their  solution  that  when  the  angle  be¬ 
tween  the  coordinate  and  laminae  natural  axes  is  degrees  in  a  Smith- 
type  plate,  the  membrane  shear  rigidity  has  no  effect  on  the  solution. 


1 


~n  Reference  -i,  the  authors  extended  the  Reissner-Stavsky  formulation  to 

sheils;  they  applied  their  solution  to  a  circular,  cylindrical,  composite 

she;:  of  Smith- type  construction  subjected  to  internal  pressure.  The 

2 

coupling  phenomenon  found  in  the  plate  analysis  was  present  also  in  that 
of  the  she! .1 .  The  authors  considered  the  validity  of  the  Donnell  approxi- 
mat ions'’  when  applied  to  composite  shells  and  concluded  that  except  for 
highly  anisotropic  composites  (E^/E^  >  1000),  any  error  introduced  into 
composite-shell  theory  above  that  introduced  into  classical  shell  theory 
is  negligible.  Whitney  and  Leissa^  reformulated  the  anisotropic  plate 
problem  in  terms  of  the  inplane  and  lateral  displacements,  and  the  stress 
resultants.  They  examined  several  bending,  stability,  and  vibration 
problems  not  previously  considered. 

The  principal  thrust  of  the  foregoing  studies  was  to  show  the  effects  of 
laminae  anisotropy  and  variations  in  laminae  orientations  on  the  behavior 
of  composite  plates  and  shells.  As  all  of  the  solutions  relied  on  the 
Kirchhoff-Love  hypothesis,  the  effects  of  interlaminar  shear  were  neglected 
entirely.  With  regard  to  interlaminar  shear  effects,  Paganob  pointed  out 
their  importance  in  discussing  the  differences  between  the  flexural  and 
extensional  moduli  of  composite  materials.  He  used  an  approximate 
method  of  Hoff  (see  Reference  6)  in  calculating  the  effect  of  shear  on 
lateral  deflections  for  the  special  case  of  bidirectional  composite  beams 
in  cylindrical  bending. 

A  more  rigorous  analysis  of  interlaminar  shear  effects  in  the  bending  of 

H 

composite  plates  was  presented  by  Whitney.  This  solution  was  based  on 
the  anisotropic  plate  equations,  including  the  effects  of  transverse  shear, 

O 

developed  by  Ambartsumyan.  Whitney  modified  the  Ambartsumyan  equations  to 
accommodate  composite  plates,  and  he  also  included  coupling  between  the  two 
transverse  shearing  strains  not  present  in  the  Ambartsumyan  formulation. 

The  solutions  obtained  by  Whitney  were  applicable  primarily  to  composite 
plates  symmetrical  about  a  median  plane;  for  unsymmetrical  constructions 
he  considered  only  a  plate  strip.  In  addition,  Whitney  s  results  were 
confined  to  the  consideration  of  the  effects  of  large  length-to-thickness 
ratios. 


2 


.  i  o 

By  extending  the  work  of  Khot/  to  include  a  first  approximation  to  the 
effects  of  inter] aminar  shear,  Taylor  and  Mayers1***  obtained  solutions  for 
axia.  compression  buckling  of  composite,  circular,  cylindrical  shells. 

In  their  analysis,  Taylor  and  Mayers  used  a  modified-Reissner  variational 
principle11  in  conjunction  with  the  constitutive  and  strain-displ acement 
relations  given  in  Reference  3  and  the  curvature  relations  analogous  to 
those  used  in  Reference  12.  Results  were  obtained  for  boron-  and  glass- 
epoxy  cylinders. 

i  110  theoretica.  development  and  analysis  presented  herein  take  an  essen¬ 
tia.,  y  different  approach  from  all  of  the  foregoing.  These  formulations 
were  based  on  classical  plate  theory  in  that  they  used  the  equations 
of  motion  in  terms  of  stress  resultants  and  the  Kirchhoff-Love  hypothesis 
to  establish  the  governing  equations.  By  contrast,  the  present,  analysis, 
suggested  by  the  sandwich-plate  developments  given  in  References  13  and  lk, 
makes  use  of  nonclassical  kinematics  and  a  variational  principle  to  estab¬ 
lish  the  governing  equations  and  associated  boundary  conditions.  The  total 
potential  energy  developed  in  a  composite  plate  during  bending  is  formulated 
in  terms  of  the  displacements,  and  the  variation  with  respect  to  the  in¬ 
dependent  displacement  functions  yields  the  governing  differential  equations 
and  boundary  conditions  in  accordance  with  the  potential  energy  principle. 

The  model  used  is  a  generalization  of  that  of  Hoff  in  Reference  13  for 
sandwich  plates.  Hoff's  model  consisted  of  a  core  enclosed  by  two  mem¬ 
brane  face  sheets.  The  core  was  considered  to  extend  to  the  median  planes 
of  the  two  face  sheets  and  to  carry  only  transverse  shear.  The  behavior 
of  two  adjacent  laminae  in  a  composite  plate  can  be  represented  by  a 
similar  model  under  the  following  assumptions: 

The  laminae  behave  as  anisotropic  membranes  with  all  of  the 
direct  stress  carrying  material  located  at  the  median  surfaces. 

The  matrix  material  between  adjacent  median  surfaces  carries 
al of  the  t  ransverse  shear. 

Hence,  the  composite  plate  is  considered  to  be  a  multilayered  sandwich 
plate  with  the  matrix  material  acting  as  cores  and  the  matrix-fiber 
combination  acting  as  the  membranes.  It  should  be  noted  that  the  matrix 
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is  considered  to  carry  two  types  of  stresses  which  do  hot  couple:  thp 
matrix  acts  in  conjunction  with  the  fibers  to  transmit  the  laminae  stresses, 
and  it  acts  alone  to  resist  the  stresses  due  to  transverse:  shear. 

1  i  .  .  i  » 

•'  I  i 

For  ease  of  analysis,  compared  with  the  methods,  for  example,  of  References 
2,  3,  and  4,  a  vector-matrix  technique  is  introduced.  The  inplane  displace- i 
ments  of  the  laminated  plate  are  described  by  two  sets  of  functions, 
and  v^  (k  =  1,2, ...,n)  ,  where  k  designates  any -one  of  the  n  liyers. 

These  sets  of  functions  are  considered  to  constitute  two  n-dimensional  vec¬ 
tors  over  the  vector  field  of  the  plate  reference  surface,  and  the  total 
potential  energy  is  formulated  in  terms  of  these  two  vectors  land  the  lat-  • 
eral  displacement  function,  w  .  Approximations  to  arty  degree  of  accuracy 

of  the  two  vectors  are  then  introduced  by  means  of  linear  transformations. 

» 

These  approximations  give  the  two  sets  of  displacement  hinbtions  in  terms 
of  smaller  sets  of  generalized  displacement  functions.  The  Euler  equations 
and  all  associated  boundary  conditions  are  then  developed  in  terms  of  the 
approximated  vectors  and  the  lateral  displacement.  Actually,  'the  Euler  '  , 
equations  corresponding  to  the  two  vector  quantities  are,  in  effect,  sets  . 
of  scalar  equations  involving  the  generalized  displacement  functions  de¬ 
scribing  the  plate  behavior.  The  theory  is  not  restricted  to  the  preserva¬ 
tion  of  plane  sections  and  can  be  extended  readily  to  include  the  definition 
of  a  separate  lateral  deflection  function  for  ehch  lamina.  This  approach 
was  adopted  in  Reference  15  in  studying  the  bending  and  buckling  of  layered 
beams.  Results  of  the  analysis  for  several  bending  and  budfling  problems 

*,  I 

of  plate  and  shell  elements  are  presented  in  the  form  of  charts.  Details  > 
of  all  theoretical  developments  appear  in  the  appendixes.  , 
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STATEMENT  OF  PROBLEM  AM)  BASIC  ASSUMPTIONS  s 

The  problem  considered  is  the  effect  of  interlaminar  shear  on  the  bending 

I 

and  buckling  of  composite,  flat  plate's  and  circularly  curved  cylindrical 
shells .  ■'  i  •  • 

i  •  i 

l  ;  ;  I 

i  , 

The  composite,  structures  ar,e  considered  io  consist  of  n  laminae  of 

'  :  '  \  i  . 

equal  thickness,  with  each  lamina  composed  of  a  matrix1  of  finite  shear’ 
transmissibility  reinforced  by  continuous  fibers  in  two  perpendicular  di- 

!  i  '  i 

rections.  ■  The  fiber,  orientation  i^ay  be  different  for  different  layers. 

The  analytical  modfel  (see  Figurb  l)  consists,:  of  n  homogeneous,  aniso¬ 
tropic  laminae.  Two  tangent -plane  displacement  .functions,  u^.  and  v.  i 
(k  =  1,2, ...,n)  are  assigned  !to  each  lamina.  In  addition,  jthete  is  a 

I  '  * 

lateral  displacement  function,  w  ,  common  to  all  n  laminae .  'The'  , 
radius  of  curvature,  R  ,  lis  considered  to  be  the  same  for  all  laminae 
and  equal  to  the1  radius  of  the  lamina tp  median  surface.  1  ■ 

,  !  ' 

•'  ,l  J  , 

Each  lamina  is  assumed  to  be  in  a  State  of  plane  stress.  In  'the  trans-  • 

I  ,  !  >  '  | 

verse  planes,  only  shearing  stresses  are  considered;  all 'other  Stresses 

'  .  i 

are  assumed  to  be  negligible.  The  transverse  shearing  stresses  are  .  , 

assumed  to  be  carried  by  the  matrix  material.  No  restriction  is  placed 

1  1  i  1 

on  the  preservation  of  plane  sections.  Finally,  the  lateral  deflections 
of  the  composite  structure  are  assumed  to  be  small  in  comparison  with  its 

I  , 

overall  thickness. i  •  ■ •  ■  *  , 

i 
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STRAIN-PIS  PLACEMENT  RELATIONS  FOR  A  LAMINA  1 

1  f  I  ]  1 

—  1  5  i  | 

For  each  lamina,  when  the,  Donnell  approximations  are  used,  the  strain-' 

displacement  relations  can  be  expressed  as 

i  '  : 
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where,  as  shown  in  Figure 

-  a  .  -  4 31 

1 ,  the 

radial  displacement 

is  positive  when 

directed  inwardly. 


TRANSVERSE  SHEAR  STRAIN-DISPLACEMENT  RELATIONS  BETWEEN  ADJACENT  IAKIIIAE 

As  derived  in  Appendix  I,  the  transverse  shear  strains  between  adjacent 
.laminae  are  given  in  terms  of  the  inplane  and  lateral  displacement  funo 
tions  of  the  k  and  k  +  1  laminae  as 

oo  _  Sw  ) 
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CONSTITUTIVE  RELATIONS  FOR  A  LAMINA 


The  relations  between  the  inplane  stresses  and  strains  are  those  which 
govern  the  behavior  of  a  homogeneous,  anisotropic  medium  in  a  state  of 
plane  stress.  In  matrix  notation,  these  relations  are 
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where  the  matrix  elements,  C.  . 
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constants  by- 


are  related  to  the  engineering  material 
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POTENTIAL  ENERGY  FORMULATION 
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}  (4) 
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(5) 


Strain  Energy 

The  strain  energy  stored  in  a  laminated,  circularly  curved  plate  consists 
of  the  strain  energy  associated  with  the  plane  state  of  stress  of  the 
laminae  as  developed  in  the  classical  manner  from  Equations  (l)  and  (3), 
and  the  strain  energy  due  to  transverse  shear  as  derived  in  Appendix  I. 
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Thus,  the  strain  energy  is  expressed  in  terms  of  the  displacements  as 


Potential  of  Applied  Loads 

Two  types  of  loadings  are  considered  as  shown  in  Figure  2:  distributed, 
lateral  surface  ;adings,  and  compressive  edge  loadings  along  the  edges 
x  =  0  and  x  =  d  .  The  potential  of  the  applied  loads  for  the  two  cases 
is  given,  respectively,  by 


L  b 


In  the  first  of  Equations  (7),  w  is  the  total  displacement  measured  from 
the  undeformed  state,  whereas  in  the  second  of  Equations  (7),  w  is  the 
additional  displacement  which  occurs  during  a  buckling  process. 

Total  Potential  Energy 

The  total  potential  energy  is  given  by  Equations  (6)  and  (7),  with  either 
N  =0  for  bending  problems  or  p  =  0  for  axial  compression  buckling 
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problems,  as 


U  +  V  = 


;)•  '“I?-?)!? 

7/ ““’■// 1^1 


dxdy 


dxdy 


dxdy 


(8) 


TOTAL  POTENTIAL  ENERGY  IN  VECTOR-MATRIX  NOTATION 

The  total  potential  energy  as  given  by  Equation  (8)  involves  2n  +  1 
displacement  functions,  u^  and  vfc  (k  =  1,2, . . . ,n)  ,  and  the  lateral 
displacement  function  w  .  From  Equation  (8),  it  can  be  seen  that  the 
set  of  displacement  functions  is  subjected  to  the  same  differ¬ 
ential  and  integral  operations.  This  consideration  also  applies  to  the 
set  of  displacement  functions  v^  .  Consequently,  it  is  advantageous 
to  consider  these  sets  of  functions  as  constituting  two  n-dimensional 
vectors,  u  and  v  ,  and  to  formulate  the  total  potential  energy  as  given 
by  Equation  (8)  in  vector-matrix  notation.  This  allows  the  subsequent 
operations  of  variation  and  integration  to  be  carried  out  with  respect 
to  the  two  displacement  vectors  u  and  v  rather  thin  the  2n  displace¬ 
ment  functions  u^  and  v^  .  With 
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i=l,2,3 

j=l,2,3 


the  total  potential  energy  can  be  written  now  as 
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REDUCTION  OF  INDEPENDENT  DISPLACEMENT  FUNCTIONS  IN  TOTAL  POTENTIAL  ENERGY 


As  shown  in  detail  in  Appendix  II,  the  set  of  displacement  functions  u^ 
(k  -  l,2,...,n)  may  be  considered  as  related  by 


\U»y)  =  F1(x,y,zk)  (l6) 

k  =  1,2, ...  ,n 

where  z,  is  the  distance  from  any  arbitrary  reference  plane  to  the 
^  th, 

center  of  the  k  lamina.  When  the  function  F^  is  approximated  by 
the  first  m  +  1  terms  of  its  power  series  expansion  in  z  as 


F1(x,y,z)  -  u"(x,y)  +  zu^(x,y)  +  ...  +  zV^x^) 


(17) 


the  displacement  functions  given  by  Equation  (l6)  become 


■  +  vr +  • 


m  * 
.  +  z,  u 
k  m 


k  =  1,2, ... ,n 


(18) 


Similarly,  the  set  of  functions  v^  can  be  approximated  by 


*  *  m  * 

v.  =  v~  +  z,  v..  +  . . .  +  z.  v 
k  0  k  1  km 


k  =  1,2, ... ,n 


(19) 


Equations  (l8)  and  (19)  can  effect  a  significant  reduction  in  the  number 
of  independent  displacement  functions  for  multilayered  plates  with  n  »  1  ; 
these  two  sets  of  equations  can  be  expressed  in  vector-matrix  notation  as, 
respectively, 


-*  ~  ->  * 
u  =  Q  •  u 


(20) 


and 


v  =  Q  •  ^ 


(21) 


14 


where 


J 


A 


(26) 


-v 


,(27) 


and 


i  =  1,2,3 
j  =  1,2,3 


(28) 


the  total 
terms  of 

U  +  V  = 


+ 


4* 


potential  energy  as  given  by  Equation  (15)  can  be  written  in 

V  *  V  * 
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where  cognizance  has  been  taken  of  the  identities 


1  '  l 
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1  * 

■  cT  • 
u  > 

Again,  for  purposes  of, the  present  analysis,  it  must  be : noted  that  in 
Equation  (29),  N  =0  when  bending  problems  are  being  considered  and 
p  =  0  'when  axial  compression  buckling  problems  are  being  considered. 

I  ‘  ' 

GOVERNING  EQUATIONS  : AND  BOUNDARY  CONDITIONS  FOR  LAMINATED  ■  FLAT  (R  -»  e°) 

- - - ’■  "  —  -  1  1  - —  - — -  —  -  -  '  • 

AND  CURVED  PLATES 

i 

Extremization  of  the  total  potential  energy  functional  given  by  Equation 

i  !  *  1  * 

(29)  with  respect  to  w  ,  u  and  v  .  results  in  the  following  equi- 

1  ■  1  , 

librium  equations  and  associated  boundary  conditions.  The  details  of  the 

,  1  1  1 

application  of  the  potential  energy  principle  are  given  in  Appendix  III. 
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Eu^er  Equations 

,  ~  _  w  '  t 

t  l  .  C22-.  !  .?-r 
1  R  R 


1 


! 


J 


(31) 


i 


I 


I  » 


I 


17 


j 


t  Cu 


«»* 

-  t  c33 


( 


2-»  *  * 
O  U  O  V 


dy  dxdy 


) 


***#■ 

‘  1  Ci2 


-,2_»  * 

o  v 


dxdy 


- 

' 1  Cl3 '  &T~  13 '  — 


dxdy 


^2-»  * 
o  v 


Sw 


c23 

by2 

G  —  B 
Sx 

+  -  A  •  U 

t 

dw 

5x 

r12.- 

8w  ^ 
by  *23 

=  0 

*2-.  * 

O  V 

+  n  .  1 

.  .2-,  *  v2_»  *  v 

p  v  (  a  u  v 

C22  ' 

by2 

t  c33  | 

\  bx2  bxby  ) 

* 

o  u 

v2_»*  * 

a  u  ^  a  U 

+  r  «  _ _ — 

C13  ' 

’  dx2  ' 

t  C12 

SxSy  23  ^ 

bv 

#*  +  ^r  • 

t 

'  l23 

bxby 

Gb  ay 

(32) 


t  bv 


t  bv 


+  —  K„3  +  —  Kpp  ® 

R  dx  23  R  5y 


(33) 


Equations  (31),  (32)  and  (33)  are  the  governing  equations,  expressed  in 
terms  of  generalized  displacements,  corresponding  to  equilibrium  in  the 
directions  of  the  displacement  w  and  the  generalized  displacement  vec¬ 
tors  u  *  and  v  *  ,  respectively.  For  application  to  bending  problems, 

U  is  set  equal  to  zero  in  Equation  (31),  whereas  for  the  eigenvalue 
problem  of  axial  compression  buckling,  p  is  set  equal  to  zero. 
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Boundary  Conditions 
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(34) 


(35) 
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(37) 


=  0  (38) 
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At  x  =  0,L 


Equations  (34)  and  (35)  require  that  either  the  transverse  shearing  force 
or  the  lateral  displacement  vanishes  at  the  plate  edges.  Equations  (36) 
and  (37)  require  that  either  the  generalized,  axial  force  vector  or  the 
generalized  displacement  vector  in  the  direction  of  the  axial  force 
vector  vanishes  at  the  plate  edges.  Finally,  Equations  (38)  and  (39) 
require  that  either  the  generalized  shearing  force  vector  or  the  gen¬ 
eralized  displacement  vector  in  the  direction  of  shearing  force  vector 
vanishes  at  the  plate  edges. 

In  viewing  the  boundary  conditions,  Equations  (34)  through  (39),  it  can 
be  seen  that,  unlike  conventional  plate  or  shell  boundary  conditions,  there 
is  no  condition  at  any  edge  corresponding  to  either  arbitrary  moments  or 
arbitrary  rotations.  The  mathematical  model  is  such  that  the  bending 
moment  at  an  edge  is  represented  by  the  couple  of  direct  stresses  acting 
through  the  n  laminae.  Thus,  the  bending  effect  appears  indirectly  in 
the  boundary  conditions  given  by  Equations  (36)  and  (38).  For  example, 
in  the  case  of  simple  support  there  should  be  no  resultant  bending 
moment  formed  by  the  laminae  membrane  stresses  at  the  edges. 

The  governing  equations  and  boundary  conditions  given  by  Equations  (31) 
through  (33)  and  (34)  through  (39),  respectively,  have  been  developed  for 
laminated,  composite,  anisotropic  plate  structures  characterized  by  a  matrix 
of  finite  transverse  shear  rigidity  reinforced  by  high-strength  fibers.  For 
application  to  the  case  of  plates  with  interlaminar  shear  effects  neglected, 
the  governing  equations  and  boundary  conditions  should  be  used  under  the  con 
dition  that  Gg  -»  E^  .  The  specialized  kinematics  of  the  problem  preclude 
the  case  Gg  -»»  unless  Equations  (31)  through  (33)  are  coupled  by 
linear  operators  into  a  single  governing  partied,  differential  equation 
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(of  fourth  order  for  plates  and  eighth  order  for  shells)  and,  con¬ 
comitantly,  the  boundary  conditions,  Equations  (3*0  through  (39), 
are  suitably  reduced. 
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METHOD  OF  SOLUTION 


For  problems  wherein  the  plate  or  shell  properties  are  orthotropic,  with 

— ) 

the  natural  and  coordinate  axes  coinciding,  the  vector  and  matrices  , 

C*  and  cl .  vanish  in  both  the  governing  equations  and  boundary  condi - 

13  23 

tions.  Consequently,  for  either  a  uniform,  compressive  edge  loading  or  a 
sinusoidally  distributed  surface  loading  acting  on  a  flat  or  curved  plate 
with  classical  simple- support  boundary  conditions,  the  displacement  func¬ 
tions 
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satisfy  automatically  the  governing  equations  and  boundary  conditions. 
Substitution  of  Equations  (40)  into  Equations  (31)  through  (33)  leads 
to  a  system  of  algebraic  equations  which  may  be  solved  for  either  bending 
(N  =  0)  or  buckling  (p  =  0)  problems,  as  appropriate. 


For  laminae  in  which  the  natural  and  coordinate  axes  do  not  coincide 
(anisotropic  media),  the  vector  and  matrices  >  ^13 >  an(^  ^23  ao 

not  vanish;  thus,  the  complexity  of  the  equilibrium  equations  precludes 
a  straightforward  solution,  even  with  the  usually  convenient  assumption 
of  simple- support  boundary  conditions.  Solutions  for  this  class  of 
problems  are  obtained  by  applying  the  direct  method  of  the  calculus  of 
variations  to  the  total  potential  energy  functional  (Rayleigh-Ritz 
procedure).  The  minimizing  sequences  are  taken  as 
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where  each  term  of  the  sequence  for  w  satisfies  the  geometric  boundary 
conditions  of  classical,  simple  support.  Details  of  the  procedure  for 
solving  the  problems  of  bending  and  buckling  of  anisotropic  flat  plates 
are  given  in  Appendix  IV. 


Equations  (40)  are  used  in  conjunction  with  Equations  (31)  through  (33) 
to  effect  solutions  for  the  bending  of  simply  supported,  square,  com¬ 
posite  plates  with  orthctropic  laminae  subjected  to  the  loading 
p  =  pQ  sin  ^  sin  ^  .  Solutions  are  obtained  corresponding  to  several 

values  of  the  parameter  ir.  in  Equation  (24),  and  the  results  are  dis¬ 
played  in  Figure  3*  Equations  41)  are  used  in  conjunction  with 
Equation  (29)  to  obtain  solutions  for  square,  flat  plates  with  aniso¬ 
tropic  laminae  subjected  to  either  a  uniformly  distributed  surface 
loading  or  a  uniform,  compressive  loading  in  the  axial  direction  (see 
Figure  2);  the  results  are  given  in  Figures  4  through  9»  Equations  (40) 
and  (31)  through  (33)  are  used  again  to  obtain  solutions  for  the  buckling 
of  infinitely  long,  flat  and  circularly  curved  plates  subjected  to  a 
uniform  compressive  loading  in  the  axial  direction  (see  Figure  2),  and 
the  results  are  presented  in  Figures  10,  11,  and  12. 
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Figure  3.  Flat  Plate  Maximum  Deflections  Corresponding  to  Various 
Approximations  to  Thickness  Distributions  of  u^  and  vk- 
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Figure  8.  Effects  of  Interlaminar  Shear  on  Axial  Compression  Buckling 
of -Flat ; Plates. 
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•’igure  10.  Effects  of  Interlaminar  Shear  on  Axial  Compression  Buckling 
of  Long,  Circularly  Curved,  Cylindrical  Plates. 


RESULTS  AND  DISCUSSION 


The  theory  developed  In  this  study  is  intended  for  application  to 
laminated  plates  and  shells.  However,  the  equations  presented  can  be 
modified  easily  and  applied  to  a  much  wider  class  of  problems.  In 
the  development,  plates  and  shells  are  considered  to  be  composed  of  a 
number  of  thin  laminae;  however,  there  is  no  requirement  that  there  be 
a  one-to-one  relationship  between  the  laminae  in  the  actual  structure 
and  those  of  the  analytical  model.  Thus,  the  complex  behavior  of,  for 
example,  conventionally  thick  plates  and  multilayered  sandwich  construc¬ 
tions  can  be  approximated  by  representing  the  structures  as  laminates. 
Individual  lamina  can  be  assigned  properties  appropriate  to  the  physical 
description  of  the  actual  structure.  The  kinematics  of  the  laminae 
provide  automatically  for  transverse  shear  deformations.  Then,  by 
introducing  independence  of  the  lateral  displacements  of  each  lamina, 
the  effects  of  both  transverse  shear  and  normal  strain  can  be  taken 
into  account.  Ordinarily,  such  modeling  would  imply  extensive,  if  not 
prohibitive,  computations.  However,  by  the  introduction  of  the  vector- 
matrix  techniques  used  in  this  study,  the  apparent  tediousness  of  analysis 
can  be  greatly  reduced. 

APPLICATION  TO  LAMINATED,  COMPOSITE  FLAT  AND  CURVED  PIATES 
Validity  of  "Plane  Sections"  Assumption 

The  degree  of  approximation  of  the  displacement  vectors  is  determined 
by  the  number  of  columns  in  the  matrix  Q  appearing  in  Equations  (20) 
and  (21)  and  given  by  Equation  (24).  In  the  latter  equation,  the  number 
of  columns  in  Q  is  indicated  as  m  +  1  ,  where  m  is  left  arbitrary. 

The  effect  of  the  matrix  Q  on  the  analysis  is  to  impose  a  constraint 
on  the  distortion  of  a  transverse  cross  section.  If  m  were  set  equal 
to  unity,  cross  sections  would  be  constrained  to  remain  plane  (but  not 
normal)  after  deformation.  Higher  values  of  m  allow  cross  sections  to 
distort  into  shapes  described  by  higher  degree  polynomials. 
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The  choice  of  a  suitable  value  of  ra  for  use  in  the  analyses  undertaken 
herein  is  based  on  the  results  presented  in  Figure  j.  There,  the  curves 
show  the  effects  of  the  shear  stiffness  parameter,  G^/E.  (  ,  on  the  de¬ 

flections  of  a  laterally  loaded,  simply  supported,  square  plate  as  given 
by  solutions  corresponding  to  various  values  of  m  .  The  ordinates  in 
Figure  3  are  the  maximum,  normalized  deflections,  with  normalization 
carried  out  with  respect  to  the  deflections  of  a  similar  plate  having 
infinite  transverse  shear  rigidity.  The  plate  properties  are  indicated 
in  the  figure;  these  properties  tend  to  maximize  the  effects  of  trans¬ 
verse  shear  in  plate-type  structures.  The  plate  loading  is  of  the  type 
jtx  .  jty 

p  =  p0  sin  —  sun  — . 

Solutions  corresponding  to  three  values  of  m  in  Equation  (24)  are 
shown  in  Figure  3.  In  addition,  the  solution  corresponding  to  the  case 
in  which  no  approximation  is  introduced  is  shown  and  labeled  m  =  n. 

The  curves  of  Figure  3  illustrate  that  only  a  small  increase  in  accuracy 
is  obtained  by  changing  from  m  =  1  to  m  =  3  ,  and  that  no  appreciable 
increase  in  accuracy  is  obtained  in  taking  m  >  3  •  In  all,  the  as¬ 
sumption  of  transverse  cross  sections  remaining  plane,  implied  in  using 
m  =  1  ,  results  in  a  good  approximation  of  the  effects  of  interlaminar 

shear  for  the  parameters  given  in  Figure  3.  Thus,  at  least  on  the  basis 
of  the  results  obtained  here,  the  assumption  of  plane  sections  remaining 
plane  but  not  normal  used  in  Reference  9  appears  to  be  valid.  The  same 
assumption  (corresponding  to  m  =  l)  is  the  basis  of  all  subsequent  cal¬ 
culations  presented  in  this  study.  The  governing  equations  and  boundary 
conditions  for  flat  plates,  Equations  (31;  through  (39),  are  expanded 
into  scalar  form  for  the  case  m  =  1  and  are  shown  to  reduce  to  the 

Reissner  plate  theory  for  homogeneous  plates  with  the  effects  of  trans- 

17 

verse  shear  included. 

Coupling  Between  Bending  and  Membrane  Stresses 

The  coupling  between  bending  and  membrane  stresses  noted  in  References  2 
through  4  is  present  in  this  analysis  in  the  form  of  coupling  between  coef¬ 
ficients  of  the  even  and  odd  powers  of  z  (see  Equation  24'.  The  terms 
involving  even  powers  of  z  correspond  to  extensional  displacements, 
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wh  jiean  terms  Lnvo.lv.!  tig  odd  powers  of  2  eorror  pot  id  to  bending  displace¬ 
ments.1^  In  the  case  of  structures  which  are  symmetric  about  a  median 
surface  and  have  the  reference  surface  (2  =  0)  coinciding  with  this 
median  surface,  the  even-powered  terms  in  Equation  (24 )  can  have  no  net 
effect  on  the  results.  For  laminates  symmetric  about  a  median  surface, 
but  having  the  reference  surface  (2  =  0)  at  any  other  location,  the 
even-  and  odd-powered  terms  couple  in  any  solution.  This  is  the  case 
with  the  problem  results  presented  in  Figure  3  where  the  reference  sur¬ 
face  has  been  taken  at  the  median  surface  of  the  extreme  lamina.  For 
nonsymmetric  laminates  (all  other  cases  considered  in  this  section), 
coupling  is  always  present  unless  a  neutral  surface  can  be  identified 
readily  and  the  reference  surface  located  thereon.  In  the  present 
vector-matrix  approach,  no  special  consideration  need  be  given  the 
question  of  coupling  between  membrane  and  bending  stresses  as  the 
operations  involved,  unlike  the  more  conventional  techniques  (see, 
for  example,  References  2,  3,  and  4),  introduce  no  additional  com¬ 
plexity  in  that  they  adjust  automatically  for  any  choice  of  reference 
surface  .location. 

Effects  of  Parameters  GjnA'll  »  E22^E11  ’  and  V12^E11  on  Ben(iin6 
of  Plates  Under  Lateral  Loading 

The  effects  of  variations  in  the  parameters  gj_2^E11  ’  E22^E11  *  and 
v12'^E11  on  maxiinum  deflections  of  a  uniformly  loaded,  square  com¬ 
posite  plate  have  been  studied,  and  the  results  are  presented  in  Figures 
4  through  6.  In  each  figure,  curves  corresponding  to  three  different 
fiber  orientations  are  shown.  Curve  (l)  in  each  case  corresponds  to  a 
construction  in  which  the  fibers  are  aligned  parallel  to  the  plate  edges 
and  change  orientation  by  90°  increments  between  adjacent  laminae. 

Curve  (2)  corresponds  to  a  construction  in  which  the  fibers  are  alter¬ 
nately  aligned  parallel  to  the  plate  edges  and  the  plate  diagonals; 
thus,  the  fibers  change  orientation  by  45°  between  adjacent  laminae. 

Curve  (3)  corresponds  to  a  construction  in  which  the  fibers  are  always 
aligned  along  the  plate  diagonals  and  change  orientation  by  90°  between 
adjacent  laminae.  The  other  properties  are  as  indicated  in  the  figures. 
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Figure  4  displays  the  effect  of  variations  in  G^2/E^  on  maximum  plate 

deflections  and  shows  that  when  the  fibers  are  oriented  along  the  plate 

diagonals  in  a  square  plate,  the  magnitude  of  G^2  has  no  effect  on  the 

deflections.  This  is  the  same  phenomenon  noted  by  Reissner  and  Stavsky 

in  considering  a  Smith -type  plate  wherein  the  angle  between  the  natural 
and  coordinate  axes  was  45°.  The  explanation  is  that  for  a  uniformly 
loaded,  square  plate,  the  directions  of  the  principal  stresses  are  par¬ 
allel  to  the  diagonals;  hence,  when  the  fibers  are  parallel  to  the 
diagonals,  the  laminae  transmit  no  shearing  stress  and  the  shear  modulus 
has  no  effect.  For  the  other  two  types  of  construction,  Figure  4  shows 
the  ratio  of  laminae  shear  to  extensional  moduli  to  have  a  considerable 
effect  on  the  plate  behavior.  It  is  evident  that  constructions  of 
type  (l)  are  particularly  poor  choices  for  composites  having  low  ratios 
of  G^2/E-q  •  For  composites  constructed  in  the  manner  of  type  (l), 

this  effect  can  diminish  the  potential  advantage  in  using  very  high 
strength  fibers.  It  is  interesting  to  note  that  at  approximately 
Gi2/Eh  =0.4  the  relative  rigidities  of  the  constructions  reverse. 

For  ratios  higher  than  0.4  the  ability  of  the  laminae  to  resist  inplane 
shear  exceeds  its  ability  to  resist  extension  and,  hence,  the  most 
favorable  construction  is  that  wherein  relatively  high  shearing  stresses 
can  be  induced. 

The  essentially  parallel  nature  of  the  curves  in  Figures  5  and  6  indicates 
that  while  the  ratios  ^22^11  and  V12^E11  of  course>  have  an  ef¬ 

fect  on  the  bending  rigidity,  this  effect  is  only  slightly  influenced  by 
the  type  of  construction.  The  differences  in  deflections  obtained  for 
the  three  types  of  construction  in  Figures  5  and  6  can  be  attributed, 
in  accordance  with  Figure  4,  to  the  effect  of  G-^/E^  =  0.2  . 

Effects  of  Interlaminar  Shear  on  the  Bending  and  Buckling  of  Flat  Hates 

The  effects  of  interlaminar  shear  on  the  bending  and  buckling  of  com¬ 
posite  flat  plates  are  illustrated  by  considering  the  behavior  of  square, 
simply  supported  plates  of  varying  aspect  ratio  subjected  to  either  a 
uniformly  distributed  surface  loading  or  a  uniform  axial  compression 
loading.  Results,  normalized  with  respect  to  the  behavior  of  similar 


plates  having  the  shear  stiffness  parameter  Gg/E-Q  -* 1  ,  are  given  in 

Figures  7  and  8  for  the  bending  and  buckling  problems,  respectively. 

Other  physical  constants  for  the  plates  analyzed  are  as  indicated  in 
the  figures. 

Representative  values  for  G  /e.^  are  0.1  (glass  epoxy)  and  0.025 
(boron-epoxy).  This  range  should  be  representative  also  of  some  of 
the  newer  plastic  composites  under  consideration  which  utilize  graphite, 
silicon  carbide,  and  beryllium  fibers  as  the  reinforcing  elements.  For 
relatively  thin  plates,  b/h  >30  ,  the  bending  and  buckling  are  sig¬ 

nificantly  influenced  by  interlaminar  shear  effects  when  the  shear 
stiffness  parameter  is  in  the  range  Gg/E^  <  0.04  .  On  the  other 

hand,  for  relatively  thicker  plates,  say,  b/h  <  15  ,  the  effects  of 

interlaminar  shear  are  significant  wnen  Gg/E^  <0.1  .  Of  importance 

is  the  fact  that  for  composite  plates  in  the  initial  portion  of  the  thin 
plate  range  and  Gg/E.^  <  0.025  ,  small  reductions  in  the  shear  stiffness 
parameter  produce  large  increases  in  plate  shear  flexibility.  Thus,  in 
the  bending  and  buckling  analyses  of  practical  thin  plates  and  plate  ele¬ 
ments  in,  say,  stiffeners  made  of  very  high-strength-fiber  composites,  such 
as  boron-epoxy,  interlaminar  shear  effects  should  not  be  overlooked. 

Design  data  for  composite  plates  with  properties  other  than  those  used 
to  obtain  the  results  presented  in  Figures  7  and  8  can  be  derived  from 
charts  similar  to  Figures  k,  5,  6,  and  9-  Charts  such  as  Figure  9  could 
be  used  to  first  determine,  for  example,  the  buckling  stress  corresponding 
to  prescribed  values  of  b/h  and  Gg/E^  and  a  reference  set  of  laminae 

parameters  (G12/eh  >  E22^E11  ’  V12^E11^  ’  Then>  adjustment  for  the 

actual  laminae  parameters  could  be  made  by  use  of  charts  similar  to 
Figures  U  through  6. 

Effects  of  Interlaminar  Shear  on  the  Buckling  of  Curved  Plates  and 
Shells  Under  Axial  Compression 

The  effects  of  interlaminar  shear  on  the  axial  compression  buckling  of 
long,  thin,  circularly  cylindrical  composite  plates  and  shells  having 
classical  simple- support  boundary  conditions  are  illustrated  by 
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considering  the  behavior  of  plates  having  specific  R/h  and  G^/E^ 
values  and  varying  cross-section  aspect  ratios.  The  results  are  pre¬ 
sented  in  Figures  10  and  11.  Figure  10  indicates  that  the  presence  of 
shallow  curvature  does  not  significantly  modify  the  effects  of  inter¬ 
laminar  shear  obtained  for  the  buckling  of  flat  plates  (r  _»oo)  .  For 
such  curvature,  the  long  plate  develops  only  rectangular  buckles  along 
its  length  (one  buckle  half -wavelength  in  the  circumferential  direction), 
the  buckle  configuration  obtained  in  the  case  of  flat  plates. 

In  Figure  11,  the  same  parameters  as  those  used  to  develop  Figure  10  are 
reused;  however,  a  longer  ordinate  range  is  shown.  The  results  now 
indicate  that  as  the  curved  plate  cross-section  aspect  ratio,  b/h  , 
approaches  the  value  at  which  the  plate  behaves  as  a  shell  (shell  be¬ 
havior  being  characterized  by  a  buckled  configuration  of  two  or  more 
half-waves  in  the  circumferential  direction),  the  effects  of  inter¬ 
laminar  shear  become  negligible.  This  result  means  simply  that  the 
range  of  thin  shell  behavior,  wherein  transverse  shear  effects  are 
quite  small,  has  been  reached.  When  thicker  shells  are  considered,  as 
in  Figure  12,  interlaminar  shear,  as  would  be  expected,  is  seen  to  have 
a  considerable  effect.  In  this  figure  the  critical  axial  compressive 
stress  has  been  normalized  with  respect  to  that  of  similar  shells  having 
a  shear  stiffness  parameter  value  GB/Ep^  -*  1  •  The  conclusion  regarding 

the  small  effect  of  interlaminar  shear  in  thin  composite  shells  corroborates 
the  results  given  in  Reference  10.  Therein,  Taylor  and  Mayers,  on  the  basis 
of  a  first -approximation  theory,  studied  the  effects  of  interlaminar  shear 

on  the  axial  compression  buckling  of  thin  boron-epoxy  and  glass-epoxy 

« 

cylindrical  shells;  they  concluded  that  interlaminar  shear  effects  could 
be  neglected  for  thin  composite  shells  using  current  high-performance 
fibers. 
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CONCLUDING  REMARKS 

A  small  deflection  theory  has  been  developed  to  assess  tjhe  effects  of 
interlaminar  shear  on  the  bending  and  buckling  of  flat  and  circularly 
curved  plates  of  composite  construction.  The  theoretical  developments, 
in  the  form  of  governing  equations  and  associated  boundary -conditions, 
are  sufficiently  general  to  encompass  anisotropic,  circularly  curved, 
cylindrical  plate  elements  of  nonsymmetric  cross  sectibn. '  Specific 
applications  have  been  made  to  the  bending  and  the  axial  compression 
buckling  of  simply  supported  flat  plates  and  to  the  axial  compression 
buckling  of  curved  plates  and  shells  having  classical,  simple  supports. 

!  I 

In  the  cases  investigated,  the  interlaminar  shear  effects  in  composite 
plates  having  moderate  cross  section  aspect  ratios,  say,  15  <  <  30  , 

become  significant  when  the  shear  stiffness  parameter r  G^/E^  ,  is 

less  than  about  0.0U  .  For  relatively  thicker  plates,  however,  the 
interlaminar  shear  effects  become  significant  when  the  shear  stiffness 
parameter  is  less  than  about  0.1  .  In  the  case  of  thin,  highly  curved 
plates  and  shells,  the  results  indicate  that  the  interlaminar  shear 
effects  are  negligible  when  the  shear  stiffness  parameter  exceeds  about 
0.01.  For  practical  composite  plates  and  shells  of,  for  example,  i 
boron-  or  glass-epoxy  construction,  effects  of  interlaminar  shear  are 
shown  to  be  adequately  described  by  a  first  approximation  theory 
as  that  used  in  Reference  10;  such  a  theory  is  given  by  this'  analysis 

when  the  exponent,  m  ,  in  the  power  series  approximation  for  distribu- 

*  •  »"* 

tions  through  the  plate  thickness  of  the  inplane  displacements  is  equal 
to  unity.  It  is  shown  also  that  when  m  =  1  ,  the  governing  equations 

for  flat  plates  reduce  to  the  theory  attributed  to  Reissner  for  homo¬ 
geneous  plates  with  the  effects  of  transverse  shear  included. 

Though  well  suited  for  composite  constructions,  the  present  theory  is. not 
limited  to  such  applications.  The  theory  could  be  used,  for  example,  to  1 
treat  the  complex  behavior  of  conventional  thick  plates  and  multilayered  , 
sandwich  const ructions  (for  example,  normal  strain  and  transverse  shear 
effects  in  heated  plates  with  irregular  thermal  gradients  through  the 
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thickness).  When  solutions  are  obtained  -with  the  aid  of  high  speed 
computers,  the  vector-matrix  development  used  in  the  present  analysis 

J  i  ’  t 

affords  more  than  a  notationai  advantage-'  as  the  necessary  operations 
(obviously  tedious  for  a  fully  coupled,  anisotropic  plate  problem)  can 
be  performed  by  the  computer  and,  hence,  they  need  not  :be  expanded. 

Further,  the  approximations  to  the  displacement  vectors  introduced 
effect  a  considerable  reduction  in  the  computer  time  required  to  obtain 
solutions.  For  example,  the  time  required  to  solve  the  governing  equa¬ 
tions  on  an  IBM  360,  computer  for  a  ( 20-layered  plat  6  and  m  equal  to 
unity  is  less  than  3  seconds,  whereas  the  time  required  to  solve 
the  equivalent  set  of  equations  with  no  approximation  introduced 
(m  =  20)  is  more  than  30  seconds. 

Although  the  present  theory  is  limited  to  problems  which  are  kinematically 

.  :  l 

and  constitutively  linear,  the  variational  development  can  be  extended  to 
include  nonlinear  strain-displacement  relations,  Since  many  fiber  com¬ 
posites  possess  essentially  linear  stresS-strain  relationships,!  the  ques¬ 
tion  of  nonlinear  constitutive  laws  does  not  appear  at  the  present  time 
to  pose  a  pressing  problem.  However,  the  effect  of  interlaminar  shear  on 

t 

the  behavior  off  composite  plates  and  shells  in  the  large  deflection  region  , 
could  be  quite  significant  and  should  be  assessed.  The  theoretical  analysis 
and  method  of  solution1  presented  herein  offers  an  attractive  approach  to 
this  problem. 


(< 


I  i 


LITERATURE  CITED 


1.  Smith,  C.  B.,  SOME  NEW  TYPES  OF  ORTHOTROPIC  PLATES  LAMINATED  OF 
ORTHOTROPIC  MATERIAL,  Journal  of  Applied  Mechanics,  Vol.  20, 

Trans.  ASME,  Vol.  75,  1953,  pp.  286-288. 

2.  Reissner,  E. ,  and  Stavsky,  Y. ,  BENDING  AND  STRETCHING  OF  CERTAIN 
TYPES  OF  HETEROGENEOUS  AEOLOTROPIC  ELASTIC  PLATES,  Journal  of 
Applied  Mechanics,  Vol.  28,  No.  3,  Trans.  ASME,  Vol.  83,  Series  E, 
September  1961,  pp.  14-02-108. 

3.  Dong,  S.  B. ,  Pister,  K.  S. ,  and  Taylor,  R.  L.  ,  ON  THE  THEORY  OF 
LAMINATED  ANISOTROPIC  SHELLS  AND  PLATES,  Journal  of  the  Aerospace 
Sciences,  Vol.  29,  No.  8,  August  1962,  pp.  969-975. 

4.  Whitney,  J.  M.,  and  Leissa,  A.  W.,  ANALYSIS  OF  HETEROGENEOUS  ANISO¬ 
TROPIC  PLATES,  Journal  of  Applied  Mechanics,  June  1969,  pp.  261-266. 

5.  Donnell,  L.  H. ,  STABILITY  OF  THIN-WALLED  TUBES  UNDER  TORSION;  NACA 
TR  479,  1933. 

6.  Pagano,  N.  J. ,  ANALYSIS  OF  THE  FLEXURE  TEST  OF  BIDIRECTIONAL 
COMPOSITES,  Journal  of  Composite  Materials,  Vol.  1,  1967, 

pp.  336-342. 

7.  Whitney,  J.  M. ,  THE  EFFECT  OF  TRANSVERSE  SHEAR  DEFORMATION  ON  THE 
BENDING  OF  LAMINATED  PLATES,  Journal  of  Composite  Materials,  Vol.  3, 
July  1969,  PP-  534-547. 

8.  Ambartsumyan,  S.  A.,  THEORY  OF  ANISOTROPIC  SHELLS,  NASA  TT  F-118, 

1964. 

9.  Khot,  N.  S.,  ON  THE  EFFECTS  OF  FIBER  ORIENTATION  AND  NONHOMOGENEITY 
ON  BUCKLING  AND  POSTBUCKLING  EQUILIBRIUM  BEHAVIOR  OF  FIBER -REINFORCED 
CYLINDRICAL  SHELLS  UNDER  UNIFORM  AXIAL  COMPRESSION,  AFFDL-TR-68-19, 
Air  Force  Flight  Dynamics  Laboratory,  Wright-Patterson  Air  Force 
Base,  Ohio,  1968. 


42 


10.  Mayers,  J.,  and  Taylor,  R.  M.,  A  FIRST  APPROXIMATION  THEORY  TO  THE 
EFFECTS  OF  TRANSVERSE  SHEAR  DEFORMATIONS  ON  THE  BUCKLING  AND 
VIBRATION  OF  FIBER-REINFORCED,  CIRCU1AR  CYLINDRICAL  SitELIS  - 
APPLICATION  TO  AXIAL  COMPRESSION  LOADING  OF  BORON-  AND  GIASS-EPOXY 
COMPOSITES,  USAAVIABS  Technical  Report  fo-t  ,  U.S.  Army  Aviation 
Materiel  laboratories.  Fort  Eustis,  Virginia,  March  1970,  AD  8716IO. 

11.  Mayers,  J.,  and  Nelson,  E. ,  MAXIMUM  STRENGTH  ANALYSIS  OF  POST BUCK LED 
RECTANGULAR  PIATES,  Presented  at  AIAA  Sixth  Aerospace  Sciences 
Meeting,  AIAA  Paper  No.  68-171,  New  York,  January  1968. 

12.  Libove,  C. ,  and  Batdorf  ,  S.  B. ,  A  GENERAL  SMALL-DEFLECTION  THEORY 
FOR  FIAT  SANDWICH  PLATES,  NACA  TN  1526,  1948. 

13.  Hcff,  N.  J. ,  BENDING  AND  BUCKLING  OF  RECTANGULAR  SANDWICH  PLATES, 
NACA  Technical  Note  2225,  November  1950. 

14.  Benson,  A.  S.,  and  Mayers,  J. ,  GENERAL  INSTABILITY  AND  FREE 
WRINKLING  OF  SANDWICH  PIATES  —  UNIFIED  THEORY  AND  APPLICATIONS, 

AIAA  Journal.  Vol.  5,  No.  4,  April  1967,  pp.  729.739. 

15.  Durlofsky,  H.,  and  Mayers,  J. ,  EFFECTS  OF  INTERLAMINAR  SHEAR  ON 
THE  BENDING  AND  BUCKLING  OF  LAMINATED  BEAMS,  USAAVLABS  Technical 
Report  70-7,  U.S.  Army  Aviation  Materiel  Laboratories,  Fort  Eustis, 
Virginia,  March  1970,  AD  871426. 

16.  Calcote,  L.  R. ,  THE  ANALYSIS  OF  LAMINATED  COMPOSITE  STRUCTURES, 

Van  Nostrand  Reinhold  Company,  1969. 

17.  Timoshenko,  S. ,  and  Woinowsky-krieger,  S. ,  THEORY  OF  PIATES  AND 
SHELLS,  2nd  ed. ,  McGraw-Hill,  1959. 

18.  Menlin,  R.  D. ,  AN  INTRODUCTION  TO  THE  MATHEMATICAL  THEORY  OF 
VIBRATIONS  OF  EIASTIC  PIATES,  U.S.  Army  Signal  Corps  Monograph, 
Department  of  the  Army  Project  3-99-11-022,  1955. 


43 


APPENDIX  I 

STRAIN  ENERGY  DUE  TO  TRANSVERSE  SHEAR 


For  a  matrix  material  with  infinite  shear  rigidity,  an  element  taken 
from  two  adjacent  laminae  (see  Figure  13a)  distorts  due  to  bending  as 
shown  in  Figure  13b;  the  relationship  between  the  inplane  and  lateral 
displacement  functions  is  simply 

t  .  >  :_v.i  (W) 

t 


With  finite  shear  rigidity  of  the  matrix  material,  the  cross  section 

(k) 

experiences  the  additional  rotation  yK  as  indicated  in  Figure  13c. 

xz 

Similar  considerations  apply  to  the  distortions- ip.,  the  y-z  plane. 
Thus,  the  relationships  between  both  the  inplane  and  lateral  displace¬ 
ment  functions  and  the  transverse  shearing  strains  become 
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The  strain  energy  due  to  transverse  shear  is 


•  t  ¥  //  m 

k=l  0  0  J 


dxdy 


Introduction  of  Equations  (43)  into  Equation  (44)  gives  the  final  form 
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Laminae 


of  the  strain  energy 

U 
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of  transverse  shear  as 
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APPENDIX  II 

REDUCTION  OF  THE  NUMBER  OF  INDEPENDENT  INPLANE  DISPLACEMENT  FUNCTIONS 


The  displacement  functions 
functions  of  the  x  and  y 


\  and  \  ,  k  =  1,2, ... 

coordinates  only;  that  is, 


n 


are  all 


u,. 


=  u. 


.(x,y) 


i 


k  = 


1,2, 


\  '  =  v,  (x,y)  (4? 

k  =  1,2, .  ..,n 

Due  to  the  continuity  of  the  deformed  structure,  there  exist  functions 
F^  =  F1(x,y,  z)  and  Fg  =  F?(x,y,z)  such  that 


F1(x,y,zR)  =  u^y) 


k  =  1,2, 


(V 


,n 


y  f  . 


vK(x,y) 


(49, 


where  z,  is  the  distance  from  a  reference  plane  to  the  center  of  the 
*bh  ^ 

k  n  lamina  (see  Figure  1).  The  functions  F1  'and  F2  can  be  ex¬ 
panded  into  a  power  series  in  z  .  These  expansions  are  written  as 


F1(x,y,z)  =  uQ(x,y)  +  z  u1(x,y)  +  ... 

+  zmu*(x,y)  +  ...  (50 

F2(x,y,z)  -  vQ(x,y)  +  z  v1(x,y)  +  ... 

+  x  vm(x,y)  +  ...  (511 


The  right-hand  sides  of  Equations  (50)  and  (5-1)  have  at  most  n  terms 
since,  corresponding  to  any  location  on  the  reference  surface,  F^  and 
F,,  must  pass  through  n  points  and,  therefore,  can  be  given  by  a  n-1 
degree  polynomial  in  z  . 
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Approximations  of  the  functions  and  F0  are  effected  by  terminating 

"th 

the  expansions  given  by  Equations  (50)  and  (51)  with  the  m°  term. 
Approximation  of  in  this  manner  and  introduction  of  the  approximation 

:fo  Equations  (Hf>)  give 


Equations  (52)  can  be  represented  in  vector-matrix  form  as 

~  — *  * 

u  =  Q  •  u 


(53) 


where 
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A  similar  development  for  v  results  in 


where 


v  =  Q  *  v 


r 

^0  • 


- 


m 

v  J 


n 

7 

'1 

n 


n 

\ 


(55) 


(56) 


(57) 


Thus,  Equations  (53)  and  (56)  represent  the  transformation  of  the  2n 
displacement  functions  u^  and  vk  ,  k  =  1,2,  .  ,.’,n  ,  into  the 
2(m  +  1)  generalized  displacement  functions  u£  and  v*  , 
k  =  1,2,  ...,m  .  The  transformation  is  quite  useful  in  analyzing 
the  behavior  of  thin,  laminated,  anisotropic  fl^t  and  curved  plates 
since,  normally,  m  «  n  1 
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APPENDIX  III 

VARIATION  OF  TOTAL  POTENTIAL  ENERGY 


To  establish  the  governing  equations  of  equilibrium  and  the  associated 

boundary  conditions,  a  stationary  value  of  the  total  potential  energy, 

Equation  (.'9'',  is  sought  with  respect  to  admissible  variations  in  the 

lateral  displacement  function  w  and  the  generalized  inplane  displace- 

ment  vectors  u  and  7  .  Although  the  variational  process  should, 

* 

in  general,  be  carried  out  simultaneously  with  respect  to  w  ,  u  and 
v  ,  the  formulation  of  the  present  problem  is  such  that  no  loss  of 
generality  in  the  derivation  is  incurred  by  carrying  out  the  variations 
separately. 


Variation  of  the  total  potential  energy  functional,  Equation  (29),  with 
respect  to  the  displacement  function  w  gives 
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(58) 


0  0  0  0  L 
Integration  by  parts  leads  to  the  following  scalar  differential  equation 
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and  associated  boundary  conditions. 


At  x  =  0,L 

!^w  Sw  ) 

Gn  t  (n-1) - B  •  u  -  N  —  >  Sw  =  0  (60) 

B  Sx  X  dx  ) 

At  y  =  0,b 

dw  _*  *T  *  ' 

Gn  t  (n-1) - B  •  v  6w  =  0  (6l) 

B  l  Sy  J 

Variation  of  the  total  potential  energy  functional,  Equation  (29),  with 


-»  * 

respect  to  the  generalized  displacement  vector  u  gives 
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where  the  identities  expressed  by  Equation  (30)  have  been  utilized. 

Integration  by  parts  of  Equation  (62)  leads  to  the  following  vector 
differential  equation  of  equilibrium  and  associated  vector  boundary 


conditions. 


Variation  of  the  total  potential  energy  functional,  Equation  (29),  with 
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respect  to  the  generalized  displacement  vector  gives 
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(66) 


where  the  identities  expressed  by  Equation  (30)  have  been  utilized. 


Integration  by  parts  leads  to  the  following  vector  differential  equation 
and  associated  vector  boundary  conditions. 
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(67) 


At  x  -  0,1, 


0  (68) 


(69) 


Again  it  is  noted  that  in  Equation  (59),  Nx  “  0  when  bending  problems 
are  considered  and  p  =  0  when  axial  compression  buckling  problems  are 
considered. 
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APPENDIX  IV 

METHOD  OF  SOLUTION  FOR  FIAT  PLATES  WITH  ANISOTROPIC  LAMINAE 


For  rectangular  plates  in  which  some  fibers  are  not  oriented  parallel 
to  the  plate  edges,  the  governing  differential  equations  do  not  possess 
a  closed-form  solution  consistent  with  the  boundary  conditions  of  simple 
support.  Consequently,  for  this  type  of  construction,  a  Rayleigh-Ritz 
procedure  is  used  in  conjunction  with  the  functional  given  by  Equation 
(29)  to  effect  an  approximate  solution. 


The  vectors  u  and  v*  and  the  scalar  function  w  are  approximated 
by 


with 


u  =  eiai  +  *  •••  +  Vk 

v>‘  .  r1@1  ... 


(70) 


where 


w  =  g171  +  g2r2  +  •••  +  ekrk 


atj  =  aj(x,y) 
=  Pj(x,y) 
rj  =  r3(x,y) 


(71) 
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0  0  ...  0  a±  a2  ...  ak  0  0  ...  0 
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(72) 
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and 


f  = 


11 


12 


lk 


'21 


'2k 


‘kk 


g  = 


Sc 


Equations  (?0)  can  be  written  as 


(76) 


(77) 


u  =  a  •  e 
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=  p  .  f 
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w  =  r 


(78) 


With  the  introduction  of  Equations  (78)  into  Equations  29)  with  R  -» 
and  cognizance  taken  of  the  identities,  Equations  (30),  the  total  po¬ 
tential  energy  becomes 
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Due  to  the  form  of  a  and  p  as  given  by  Equations  (72)  and  (73),  the 
matrices  <Jk  given  in  Equations  (80)  have,  in  turn,  a  special  form. 
These  matrices  are  equivalent  to  partitioned  matrices  in  which  the 


number  of  submatrices  equals  the  number  of  elements  in  C. 


or  A 


* 

and  H  ,  respectively. 


o?  = 


and  7*  as  given  by  Equation  (7*0,  the  submatrices  are  formed  by  multi- 
plying  their  corresponding  elements  in  C.  .  by  the  matrix  product  in- 
volving  the  vectors  C?  ,  ,  and  y  or  their  derivatives.  For  example, 
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2x2  matrix,  5^  is  formed  from  four  sub- 


with  Cu  assumed  to  be  a  A  ^ 

matrices  given  by  the  respective  element  of  Cn  times  the  ^atnx 


-»  #  T 

“x  •  “  x 


# 

J 


Upon  performing  the  integrations 

L  b 


?i  =  ff  k+*6  +  2ltk  + 

r\  r\  ' 


o  o 

L  b 


T?  *12  )  ** 
t  * 


i 

I 


'■  -  //  (® 


2®u  +  7 


*13  ) 


f 

dxdy 


0  0 
L  b 


-3  JJ  lV  *3  ♦*>)"* 

0  0 


0  0  ' 

L  b 

=  //  (»8  +  ^lo)4’** 


2G. 


B 


%  -  ?  If  * 


L  b 

// 

0  0 

L  b 

dxdy 


0  0 


\(83) 


sT  =  ’)**  J 

0  0  '  ! 

the  total  potential  energy  given  by  Equation  (79)  becomes 
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Extremization  of  (U  +  V)  with  respect  to  ,ef  ,•  ?  and  gf  leads  to 
■the  following  yector  equations;  '  , 
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The  characteristic  functions  y ^  and  0^  (j  =  1,2, 3, 4)  given 

by  Equations  (74),  (8l)  and.  (82),  respectively,  are,  now  taken  as  i 
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In  view  of  Equations  (80)  arid  (80),  the  integration  of  the  matrix  in¬ 
tegrals  Sj  (j  =  1,2, ..., 7)  ^  Equations  (83),  cdn  be  effected  with 
respect  to  the  functions  CTj  ,  0j  ,  and  7^  .  Subsequent  substitution 

of  the  integrated  'S.  functions  (j  =1,2,  ...,7)  in^o  the1  Euler 

tJ  l  |  1  1 

equations,  Equatibns  (05)  through  (87),  leads  to  a  set  of  algebraic' 

I  , 

equations.  The  solution  of  these  algebraic  equations  yields  the 
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desired  vectors  e  ,  f  ,'  and  g  .  1 
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APPENDIX  V 

EXPANSION  OF  GOVERNING  EQUATIONS  AND  BOUNDARY  CONDITIONS 


CORRESPONDING  TO  A  FIRST  APPROXIMATION  (m  =  l)  FOR 
FIAT  PLATES  (R  -»«);  REDUCTION  TO  REISSNER  PLATE  THEORY 


With  the  definitions 
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and  m  =  1  ,  the  governing  differential  equations,  Equations  (31)  through 

(33),  can  be  expanded  to  give  the  following  five  governing  scalar  equations 


for  the  bending  (N  =  0)  of  flat  plates  (R  -**>)  . 
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The  boundary  conditions  associated  with  Equations  (90)  through  (94)  are 
obtained  by  expanding  Equations  (34)  through  (39)  for  m  =  1  and  R  -»<» 
These  expansions  are  straightforward  and  result  in  the  following  scalar 
equations. 
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For  plates  symmetric  about  their  reference  planes,  the  material  constants 

#  * 

vanish.  Consequently,  the  membrane  displacements,  uQ  and  vQ  , 
and  the  bending  displacements,  ,  v^  ,  and  w  ,  uncouple  and  only 
Equations  (92),  (93)  and  (9^)  are  required  to  describe  the  plate  behavior. 
In  addition,  for  homogeneous  plates  (n  -» «  ,  t  -»  0)  ,  the  material 

constants  become 


Thus,  for  homogeneous  plates,  Equations  (92)  through  (91*)  reduce 


to 

o  (9fi) 


(100) 


Differentiation  of  Equations  (96)  and  (99)  with  respect  to  x  and  y  , 
respectively,  and  summation  of  the  two  resulting  equations  ..  eal 
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Then,  substitution  of  Equations  (102)  into  Equation  (101)  results  in 


DN7  v  =  - 


6(i-v) 


s/~  p  +  p 


(103: 


Equation  (103)  is  of  the  same  form  as  the  governing  equation  in  Reissr.or 

plate  theory.1^  The  only  difference  between  Equation  (lOj'  and  the 
governing  equation  of  the  Reissner  plate  theory  appears  in  the  coefficient 
of  the  first  term  on  the  right-hand  side  of  the  equation.  This  coef¬ 
ficient  is  slightly  different  in  the  Reissner  theory  because  the  trans¬ 
verse  shearing  stresses  are  assumed  therein- to  vary  parabolica  iy  through 
the  cross  section.  For  m  =  1  in  the  present  analysis,  the  transverse 
shearing  stresses  are  constant  through  the  thickness  of  the  plate. 
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